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Abstract — In this paper a general methodology is introduced for the determination of potential prototype curves used for the 
drawing of prehistoric wall-paintings. The approach includes a) preprocessing of the wall-paintings contours to properly partition 
them, according to their curvature, b) choice of prototype curves families, c) analysis and optimization in 4-manifold for a first 
estimation of the form of these prototypes, d) clustering of the contour parts and the prototypes, to determine a minimal number of 
potential guides, e) further optimization in 4-manifold, applied to each cluster separately in order to determine the exact functional 
form of the potential guides, together with the corresponding drawn contour parts. The introduced methodology simultaneously 
deals with two problems: a) the arbitrariness in data-points orientation and b) the determination of one proper form for a prototype 
curve that optimally fits the corresponding contour data. Arbitrariness in orientation has been dealt with a novel curvature based 
error, while the proper forms of curve prototypes have been exhaustively determined by embedding curvature deformations of 
the prototypes into 4-manifolds. Application of this methodology to celebrated wall-paintings excavated at Tyrins, Greece and the 
Greek island of Thera, manifests it is highly probable that these wall-paintings had been drawn by means of geometric guides 
that correspond to linear spirals and hyperbolae. These geometric forms fit the drawings' lines with an exceptionally low average 
error, less than 0.39mm. Hence, the approach suggests the existence of accurate realizations of complicated geometric entities, 
more than 1000 years before their axiomatic formulation in Classical Ages. 

Index Terms — rotation and translation invariant curve fitting, pattern recognition in paintings, optimization in differentiable 
manifolds, geometric guides in prehistoric wall paintings, minimal parameters set for curve description, fitting prototype curves to 
drawn borders. 

> 



O 

(N 

O 

On 
> 

u 

^ 1 Introduction 

In this paper a methodology of general applicability 
is presented for answering the question if an artist 
\Q used a number of archetypes to draw a painting or 
CN if he drew it by free hand. In fact, the contour line 
Q parts of the drawn objects that potentially correspond 
^—H to archetypes are initially spotted. Subsequently, the 
C^l exact form of these archetypes and their appearance 
throughout the painting is determined. The method 
^ has been applied to celebrated wall paintings, 'Tady 
of Mycenae'', of the 13th century B.C. excavated at 
Mycenae, Greece and ''Naked Boys" drawn c. 1650 
B.C. excavated at Akrotiri, Thera, Greece. 



• D. Arabadjis, P.Rousopoulos, Constantin Papapodysseus and Michalis 
Exarhos are with the School of Electrical and Computer Engineering, 
National Technical University of Athens, Iroon Polytechniou 9, 15773, 
Athens, Greece. E-mail: cpapaod@cs.ntua.gr 

• M. Panagopoulos is with Department of Audio & Visual Arts, Ionian 
University, Corfu, Greece. 

• L. Papazoglou-Manioudaki is with National Archaeological Museum 
of Greece, Patision 44, Athens 

© 2011 IEEE. Personal use of this material is permitted. Permission 
from IEEE must be obtained for all other uses, in any current or future 
media, including reprinting/republishing this material for advertising 
or promotional purposes, creating new collective works, for resale or 
redistribution to servers or lists, or reuse of any copyrighted component 
of this work in other works. 



1.1 Previous Works on Fitting Data to Curve Pro- 
totypes 

In the bibliography there are various approaches to 
the problem of optimally fitting curves to 2D-data, 
either using the explicit or the implicit form of the 
model curves. Namely, in |1|, the optimal fit of im- 
plicit plane curves to data points is treated via a 
square distance minimization procedure. This mini- 
mization is performed over translation, rotation and 
implicit function parameters, in order to deal with 
rigid body motion of the model curve and with its 
shape variance. Previously, in |2|, authors had em- 
ployed orthogonal distances and the related tangential 
quantities between data and model curve points, in 
order to optimally fit conic sections to the given data 
points, via a non-linear regression algorithm. A lower 
bound for the least squares non-linear regression, 
between an implicit model curve and the data points 
of a contour, is given in |3|, thus offering a measure 
for the optimality of the orthogonal distance curve 
fitting methods. In |4|, authors analytically determine 
the least-squares optimal translation and rotation of 
an explicit curve model, in order to fit a given set 
of data points. Then the primary parameters of the 
model curve are obtained via a 2D iterative region 
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descending process. The alignment of an implicit 
curve model to given contour points is also dealt in 
||5j, in a manner that the fitting process is invariant 
under affine transformations. Namely, the authors 
use affine invariant Fourier descriptors for an area- 
parametrization of the model curve and, by matrix 
annihilation, they determine the harmonic implicit 
form for the model curve descriptive equation. In 
1^1 authors exploit point-wise distribution of shape 
templates in order to define shape descriptors and 
develop the corresponding distance measure. This dis- 
tance measure offers exact point-wise correspondence 
between two different shapes, allowing for optimal 
matching transformation of one shape to the other. In 
||7J authors deal with the problem of fitting explicitly 
described B-Spline model curves to data-points lying 
on a 2D-surface. The B-Spline is restricted so as to lie 
on the same surface with the data-points by treating 
the surface as a differentiable 2-manifold. Then, the 
fitting procedure is performed so as the on manifold 
geodesic distance between data and model curve 
points is minimized. Analogous technical elements for 
modifying planar objects so as to lie on a manifold 
are developed in |8|, where the authors re-evaluate 
the ''Beltrami short time kernel'' in the case that 
the image to be filtered lies on a manifold. In |9| 
authors deal with the problem of fitting algebraic 
curves to data points, via their implicit representation, 
by constraining their functional form so as to have 
closed zero level sets, thus making the curve deter- 
mination more robust. The implicit functional form 
of the model curves (algebraic) is adopted in |10| 
to define itself a matching error via the polynomial 
function level sets. In order to deal with arbitrary 
orientation, authors adopted invariant computation 
techniques of polynomial coefficients. Finally, in |11| 
the representation of an object as a unification of sub- 
objects, that can be analytically described by implicit 
curve or surface forms, is used, in order to segment a 
complex visual form into "meaningful" sub-shapes. 

1.2 The Present Problem and the Proposed Ap- 
proach 

In the present work, the demand was the determina- 
tion of an unambiguous relative placement between 
a prototype implicit curve and the contour data, in 
the sense that attribution of various contours data to 
a prototype curve is orientation invariant. Namely, 
given an implicit form of a prototype curve, the 
integral of the least distances between a curve part 
and each contour data points is computed immedi- 
ately, by means of the flat curvature of the implicit 
functional form of the prototype curve, with no need 
of intermediate optimal placement of the data points 
along the prototype curve points. Additionally, for the 
present application, one must spot a unique prototype 
curve that gave rise to a considerable number of 



different contour realizations, where each such real- 
ization corresponds to a different part of the prototype 
curve. So, the determination of the model curve that 
optimally fits a set of partial realizations of its planar 
graph is a problem of increased degrees of freedom, 
compared with the widely treated problems of curve 
fitting. Namely, the problem of optimally fitting the 
model curve to the data runs over the space of the free 
parameters of the curve model and, simultaneously, 
over all possible line segments of it and all possible 
data points subsets that could have been generated 
by the same prototype. There are various approaches 
that deal with the determination of the parameters 
that optimally fit a model curve to a given set of 
points in a statistically efficient manner. In |12|, there 
is a compilation of the existing implicit curve fitting 
techniques, together with statistical analysis of their 
efficiency. Concerning the works referred and outlined 
in Sect, [l.!} they are mainly split into 2 categories: 
a)those who treat the problem of fitting curve models 
to data points and b) those who treat the problem of 
shape matching. Among these publication there are 
approaches that treat the problem of curve fitting (e.g. 
Ul/ I5l |4|, |10|), or shape matching (e.g. |^|) with 
an orientation invariant approach, but none of these 
techniques dealt with the problem of simultaneous 
fitting the same implicit curve model to different 
sets of data points. Moreover, in most curve-fitting 
approaches, determination of the optimal values for 
the free parameters of the model curve is based on 
iterative error minimization algorithms like Gauss- 
Newton, Gradient Descent, etc. that do not necessarily 
converge to a global error minimum. Concerning 
shape matching methods, although they offer unsu- 
pervised shape similarity measures, minimization of 
these measures deforms the model shape without 
constraining any of its geometrical characteristics. In 
order to deal with the demands of unambiguously 
optimal fit between the determined model curve and 
all contour data attributed to it, we had to develop 
a unified optimization method of the least possible 
dimensionality that exhaustively acts on the set of all 
possible prototypes of the same parametrized curves 
family. The fitting process introduced here treats each 
class of prototype curves as a 4-manifold and per- 
forms a curvature driven exhaustive error minimiza- 
tion between the class of prototype curves and the 
available data points. By applying this methodology 
two times, first with respect to the drawn contours 
and next with respect to the prototypes (Sects. |43)|4.6| 
respectively), we determine a unique prototype curve 
that optimally fits all data points attributed to it. 
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2 Fundamental Notions and Defini- 
tions Concerning the Method of Wall- 
Paintings Drawing 

A first crucial hypothesis for the method the artist(s) 
might employed for drawing the wall-paintings c. 3500 
years ago 

'The fresco technique asked for highly fast and pre- 
cise execution of the drawing process. At the same 
time, the stability of the contour line of the figures 
depicted on the painting had always been a primary 
goal of the artists throughout human history Thus, 
a plausible assumption about the method of drawing 
of these prehistoric wall paintings is that the artists 
employed guides or instruments to support the draw- 
ing object/' The quality of drawing of the paintings 
considered in this paper, suggests an additional con- 
tent to this assumption: ''the artists took special care 
and paid attention to ensure continuity of the drawn 
contour line and, wherever possible, of its tangent, at 
the points where a change of guide occurred." 

Definition 1 (The object part): If the aforementioned 
assumption is correct, then there will be subsets of 
the contour line of the various figures appearing in 
the wall-painting, which they had been drawn by a 
continuous stroke of the paintbrush. These subsets are 
called object parts. 

Definition 2 (The object): One can define the object 
to be a subset of the wall painting border, which 
represents a thematic unit, is smooth and its beginning 
and end points are discontinuities of the border line 
or of its tangent. Evidently an object is a concatenation 
of object parts, since it has been drawn by a usually 
small number of continuous strokes of the paintbrush. 

2.1 A rigorous determination of the object part 

Consider ideally that, in the wall painting, all ob- 
jects' contour lines are continuous, described by the 
piecewise twice differentiable mono-parametric vector 
equation f{t) = {x{t),y{t)). Then, an object part is 
spotted by determining a contiguous subset of the 
object curve whose beginning and end points are one 
of the following: l)the beginning or the end of the 
object, 2)a point where the curvature is discontinuous. 
Analogous pictorial feature definitions can be found 
in fl3\ , where authors use curvature extremes and 
sign changes in order to define strokes, corners, end- 
points and junctions. 

The notion of an object part is crucial, since if 
the artist(s) indeed used stencils, then, each time he 
placed the stencil on the wall and drew a line, we 
consider he created an object part. 

If Lo is the length of an arbitrary object then we 
select an appropriate small percentage of Lq, say L^. 
Subsequently, in order to spot points where there 
is change of stencil, we proceed as follows: Let an 
arbitrary point on object contour, say the i-th, where 
Ls ^ ^ ^ Lo — Ls, then we determine the 3rd 



degree polynomials ff (xf (5), (5)) of object length 
s e [0,Ls] that best fit object points p{i — Ls)...p{i) 
and the polynomial of 5th degree rf (xf (5), (5)) of 
object length s e [O,^^] that best fits the points 
p{i)...p{i^ Ls) and satisfies the tangent continuity on 
V{t), namely J^ff (i-j) = J^rf (0) = (±<i,y<i). 

In order to decide if the i-th pixel is an ending 
point of an object part we check for discontinuity of 
the object contour curvature on it. Hence, the two 
approximations of the contour curvature from the left 
^^{Ls) = fl{Ls) and from the right c+(0) = rf (0) 
should manifest an abrupt difference. Equivalently, 
we demand \c^{Ls) — c^(0)| > SO, where SO a very 
small angular threshold. If this demand is met then 
the i-th point p{i) is the end of the current poten- 
tial object part. The beginning of the next potential 
object part is the first point after i-th, say the j-th, 
where the demand \c~{Ls) — c^(0)| < SO is met. By 
application of this method to all available objects of 
the wall-painting figures, we generate an ensemble of 
potential object parts with contour point vectors r^, 
p = Ij ...^Mop. We once more stress that we consider 
that each such object part was drawn by one paint- 
brush stroke. We also assume that if we achieve in 
corresponding these object parts into a small number 
of prototypes, then we will support the hypothesis 
that these "one-stroke contour lines were made by 
means of a stencil or another equivalent instrument 
that guided the brush." If we further adopt the plau- 
sible assumption that each such stencil had a specific 
geometric form, namely in modern mathematics a 
specific functional type, then we may deduce that 
the derivatives of all orders of f{t) remain continuous 
along an object part. On the contrary, at the points of 
change of stencil, or of replacement of the stencil, it 
is highly probable that the curvature will essentially 
change, even in the case where the artist(s) achieve 
continuity of f{t) and f{t) to ensure a smooth and 
good aesthetic result. 

3 A Brief Description of the Em- 
ployed Methodology for Determining 
Possible Geometric Guides 

After application of the aforementioned methodology, 
we have determined an ensemble of object parts, i.e. 
an ensemble of parts of the painted figures' borders 
with small curvature fluctuations along each such 
part. Now, the question arises if there is a small 
number of prototype curves that optimally fit these 
object parts and if yes, how to determine the exact 
functional form of these prototypes. In order to deal 
with this problem we have applied a novel method- 
ology outlined in the following steps. 
Step 1 - Choice of a reasonable criterion that quantifies 
how well an object part fits a prototype curve part 
Suppose the prototype curve S in that is described 
by the implicit functional form fs{x^y) = and let 
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TABLE 1 
Basic Symbols 



Symbol Meaning 



fsiXi = implicit description of a model curve 

X (x,y) Cartesian coordinates 

a {ci,b) primary parameters of the model curve 

cs(x, oi) the flat curvature function of the S model curve cs = • nZ^'^^ii 

Xs{^, the model curve, parametrized via curve length ^, which is resulted by a fixed pair a of primary parameter values 

(^s((3, x) the primary parameters isocontour, parametrized via curve length (3, as it is resulted by 

a fixed pair x of Cartesian planar coordinates 

the curve length while moving on xy-plane in a direction normal to xs 
Sa the curve length while moving on a6-plane in a direction normal to q;^ 

the internal product between the directions normal to q;^ and xs 
{T^,T^) (\\dxfs\\,\\dcfs\\) 

{drj^ (i7, (irj, c^7^) The local frame defined by the level-sets of fs'- dij is the curve length differential along Vfs and 

the remaining curve length differentials evolve along the curves selected to define the tangent space of fs- 
dj evolves along (xs^c^s) direction, dr^ evolves in a direction normal to covariance of x and a 
and dj^ normal to covariance of x and /3 

(n,l) the unit directional vectors normal and tangent to xs respectively 

(tT'S^Is) the unit directional vector spaces normal and tangent to the level-sets of fs respectively 

^s(x^ ^) the Euclidean norm of fts evolution on the level-sets of fs (i.e. of the level-sets' curvature) 



n be the unit vectors normal to it. Moreover, let an 
arbitrary point M on an object part and its vectorial 
distance Sm from S. Suppose, for a moment, that the 
object part is a continuous curve. Then, evidently, a 
measure for the level of fitting the considered object 
part to 5* is the integral of Sm = Sm along a proper 
part of S. 

Step 2 - Alternative version of the aforementioned 
criterion 

Because Sm is calculated along n^, by applying Stokes 
theorem on domain bounded by the object part and 
the curve part of S, we obtain. 



2 Sm • ndl = / SMdl = / (V • n 
Js Jan Jn 



)dn (1) 



But, actually, V • n = c^, where cs is the flat curvature 
of the implicit functional form fs of the curve model 
S. Consequently, minimization of the integral of Sm 
on a proper part of S can be obtained by minimization 
of the corresponding area integral of cs- In Sect. 4.1 
and in appendix A, we show that minimization of the 
integral of the flat curvature over Q can be obtained 
by minimizing the integral of cs on the boundary of 

a on 

Step 3 - Primary Parameters definition 
In general, the functional form of fs depends on 
a significant number of parameters. For example, 
conic sections implicit functional form depends on 
5 parameters in general. In Sect. |4.2[ it is shown 
that, independently of the number of parameters and 
the exact functional form of fsr two parameters are 
sufficient to describe all possible deformations of the 
flat curvature of fs on the (x, y)-p\ane. We will call 
these parameters as the ''primary parameters'' of fs- 
Step 4 - Extending the fitting process in 
Assume a given object part, OP and a parametrized 



class of prototypes S, implicitly described by the 
equation fs{x,y,a,b) = 0. We look for the determi- 
nation of the proper primary parameters values (a, b) 
and the proper relative placement between OP and 
S, so as the integral of the point by point distances 
(see Step 1) between these curves are minimized. 
Hence, we extend the object part - prototype fitting 
analysis from 3?^ to 3?^, where two of the coordinates 
of refer to the {x^y) coordinates of the curves 
in 3?^, while the other two coordinates refer to the 
primary parameters (a, b) of fs- At the beginning, this 
specific subset of 3?^ is assumed to be the Cartesian 
product of the (x, ?/)-space with the (a, 6)-space. But, 
in the following, given that we need to calculate 
distances between OP and all possible deformed ver- 
sions of S, we treat the space of the level sets of 
fs as a differentiable manifold. Via the curvature of 
the tangent space of this manifold, we look for the 
curve S which optimally fits OP along the normal 
vectors fis as it is described in Sects, |4.3| 



4.4 The 



aforementioned approach and the resulting criterion 

form an efficient procedure that simultaneously offers 

the primary parameters of the prototype curve and its 

optimal relative position with OP. 

Step 5 - Determining the different types of potential 

stencil-guides 

After application of the fitting process outlined in Step 
4, each object part corresponds to a prototype curve 
part of a specific implicit functional form fs with 
primary parameters (a, 6). Then, we apply standard 
clustering techniques in order to group the estimated 
primary parameters of the same class of prototype 
curves into clusters of neighboring parameters values. 
We assume that each such cluster corresponds to a 
specific stencil-guide. In other words, we assume that 
all object parts that correspond to the same cluster 
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of primary parameters have been generated by one 
stencil - guide, the optimal cluster representative. 
Step 6 - Determination of the optimal cluster represen- 
tative, i.e. of the potential guides parameters 
Consider a cluster and the set of object parts that 
correspond to it, OP^. Then, we redefine the sub- 
domain between the level sets of the prototype curves 
S and the image of OP^ on the -manifold, as it is 
determined via the estimated (a, 6) G Q^. Application 
of the resulting criterion, (see Sect. |4.6[ ), offers a unique 
prototype curve with specific primary parameters, to 
which all OPg are optimally fitted. 

We stress that the aforementioned procedure ac- 
counts for all possible prototype curve versions, de- 
scribed by the functional form fs and acts on the 
manifold they define independently of rotation and 
translation transforms. So, the fitting process is inde- 
pendent of the initial orientation of the object parts 
on the xy -plane. Moreover, by computing the fitting 
error as an integral on the manifold of a class of 
prototype curves and by performing the error mini- 
mization exhaustively, the developed procedure offers 
very consistent fitting results, even in the case that the 
considered drawings suffer serious wear. 

4 DiFFERENTIABLE MANIFOLD ANALYSIS 

FOR Fitting the Outlines of a Painting 
TO Prototype Curves 

4.1 Determination of the most probable stencil 
part that created an object part 

Let an object part consist of Np pixels with point 
vectors rp{i)r i = l...A/p. Then the contour length Si 
that corresponds to each point vector rp{i) is approx- 
imated via Si = Yl^jZ^i \\rp{k + 1) - rp{k)\\. So, we can 
express the object part contour with curve parameter 
the length s, namely rp{si). 

If the painter had used a part of a specific prototype 
curve to draw this object part, we expect that the 
drawn contour and the corresponding part of the 
prototype curve should have close curvature values 
along them. So let a prototype curve (stencil) S given 
by the equation fs{x^y) = 0. Curvature at length s 
of this stencil cs{s) = xs{s)ys{s) - ys{s)xs{s), is also 
obtained via the relation 



c{xs{s),ys{s)) = V ■ [j^) 



ixsis),ys(s)) 



In order to compute the minimum distance of 
an object part from a properly prototype curve we 
adopt the flat version of the curvature c{x,y) = 
and for any point {xp^yp) of the 



^ ll|V/s|| 



object part we compute the value c{xp,yp). Assume, 
for a moment, that we properly place the object part 
around the adopted prototype curve so as to minimize 

between 



V/s 



sis),ysis)) 



the distances 6{s) = 5{s) 

them. For the computation of the integral of these 
minimum distances along the prototype curve we 



define the domain ft between the prototype curve and 
the object part. Then, using Stokes' Theorem |14| for 
the integral of c(x, y) on ft, we obtain 



Ii{so,S,p) = J \c{x,y)\dn = 



V/s 



l|V/s|| 



S{s) 



ds 



nSo-\-Lp 

h{so,S,p)=2 / 6{s)ds 

J So 



(2) 



(3) 



We also use Stokes' Theorem to obtain a boundary 
form for the integral of c^{x^y) on ft, as described 
below: 



hiso.S.p) = J J^c{x,y)^dn 

= n^v-{c^)dn 



l|V/5| 



(4) 



Since n = is the unit normal vector of fs{x,y) 

isocontours, the area differential of ft can be written 



dft = dSfidsl, where l^n and 
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Vc{x,y) 



l|V/; 



1. So, 



(5) 



Then, by applying Stokes' Theorem in (|4| we obtain 

nSQ-\-Lp 

h{sQ,S,p)= I \c{fp{s- SQ))-c{fs{s))\ds (6) 



In order to obtain the most possible part of sten- 
cil 5* used for drawing the considered object part, 
we should determine the proper sq that minimizes 
Ii{sq,S,p). But as it is shown in Appendix |A] min- 
imization of 12(50, 5', p) coincides with minimization 
of Ii{sq,S,p). Thus, the part of stencil S that had, 
most probably, been used for drawing the consid- 
ered object part consist of points 7% = {xs{s),ys{s)), 
s e [sq.Sq Lp], where 



sl = argminS^ \c{fp{si)) - cs{si + so)\ (7) 



4.2 Definition of the primary parameters for a fam- 
ily of prototype curves 

The analysis of previous Sect. 



4.1 



has offered a func- 
tional representation of the optimal fitting between 
an object part and a given prototype, S, described 
by the locus r G 3?^ : /s'(r) = 0. But, the prototype 
curve, S, can be described as a realization out of all 
the curves that belong to the parametrized family 
fs{x,y\A) = 0; where A = [ai, aAr^] is the vec- 
tor of the free parameters for this class of curves. 
The parameters that form vector A describe spatial 
transformations (rotation and translation) and shape 
(curvature) variances between the curves of the corre- 
sponding family. Since, analysis of Sect. |4.1| makes use 
of the flat curvature functional over the descriptive 
equation of the prototype curve, we should determine 
the number of parameters which are sufficient to 
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/Select an Object Part 
Contour rp(s) 



Select a Prototype 
Curves' family 5witti 
primary parameters 

a={ab) 



Find proper starting 
points CO*(p) on S 
ttiat satisfy equation 
(25) 



Among all fitting error 
stationary points Co*(p) find 
the proper pair {/o*,ac*){p) 
ttiat satisfies equation (23) and 
calculate the corresponding 
minimal error £?7-(5,r;;.) 



Select the prototype 
cun/es family offering 
the minimal Err(S,rp) 
and attribute Object Part 
Vp to this family 



^re Object Pa^ts^ 
exhausted 



Using all repositioned rp(s) find the 
proper point ^ that satisfies 
equation (31} of fitting error 
stationary points and obtain the 
optimal primary parameters of 
prototype curve of Qs^ by 
numerically solving equation 




For each object part rp{s) find 
the proper repositioning CO*{p) 

that aligns rp{s) so as to 
minimize error (30), which is a 

measure of fitting a single 
prototype to all elements of Qs"^ 



Select a group Qs'^ and 
/ consider all Object Parts that 

belong to it and the 
corresponding prototype curve; 
estimates {/o*,otc*){p) 




Cluster all primary parameters 
estimates {ac*)(p) into 
groups Qs'^ of neighboring 
parametrs values 



Fig. 1 . Block Diagram outline of the developed methodology 



describe all curvature variances between the curves 
of the same family. 

Namely, we should exploit the derivative 
da^csix^A), where x = {x,y)^, ^ = {-y,x)^±Xf 
and cs{x^A) is the flat curvature that corresponds 
to the class of curves with descriptive equation 
fs{x\A.) = 0. For the subsequent analysis we should 
recall the definition of the unit directional vectors 

^j|^±n(x) employed in 
Sect. ^ ^ ^^ 
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where V^^ = (dx.dy) , V^^= {-dy.dx) • So, 
the flat curvature is related with l{x) arid n{x) via 
the expressions 



_9_ 
dx 



)T rj-i 
i=-csix,A)i (l^f) i=csix,A)h (8) 



and the partial derivatives of the directional vectors 
with respect to parameter a„ read 



|Vv/s 



I 



Subsequently, using expression (|8| for the curvature 
and (|9jl for da j and da„n, the a„-derivative of cs reads 

da„csix,A) = da 



[da J 

for 



iS) and {9) imply that 9a J = (cstj 



n] L, smce 



= 0. Using equation (9| for da^n and next (|8| 
/ we obtain 



IIVxMI 



Cs \n 



WxfsW 



(10) 



Next, by expanding the derivative, V^^, of the function 
fs^ = dajsr one obtains V^/^- = a^/^^n + dits^l 
and subsequently (IlO| results 



IVv/sll 



(11) 



Using equation ( pT| we obtain the differential of 

cs{x-,A) in the parameters space d^cs = {d^cs)dA = 



^^11 rV^{l V^d^fsdA), since / V^^ commutes 
with addition. But the maximal number of linearly 
independent vectors S/^d^fsdA is 2. Equivalently the 
minimal number of parameters sufficient to describe 
all curvature variances among the curves of the same 
family is 2. 

In order extract the primary parameters (a, b) of a 
curves family parametrized via a parameters vector 
A of arbitrary size, we exploit the influence of affine 
mapping to the flat curvature function. Namely, if x 

{A^R) XT + d, where R : R^R = /, Ax ^ ^ 



b 



csixT.A) = {ab)^cs{x.A) 



(12) 



Next, since a, b are linearly independent, 
we suppose that any parameters (i.e. non- 
spatial) variation of cs can be obtained 
by variations of (a, 6), da = (da.db). So, 



_{abr 



and since du Wd^fsl 



lAxV^/sll 



da 



dA 



cs{x,A) 
'dii\\d^fs\ 



Vi In 



|AxVx/s| 



da (13) 



where = and On = Fv^(Fv^). 

Hence, given a parametrized implicit curve 
fsiXiA) = 0, its primary parameters a = {a,b) can 
be recovered by the mapping {x,y) ^ {x/a,y/b) 
and then any differential variation of the primary 
parameters evolves free parameters vector A as 
equation (Il3| describes. 
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4.3 Analysis for optimally fitting an object part to 
a class of prototype curves 

In Sects. |4A] and 14.21 it has been shown that the action 
of the flat curvature functional over the implicit form 
of a class of curves fs{x,y,a,b) = is sufficient to 
optimally align an object part along the curve defined 
by a fixed parametrization (a, 6) = A. But, still, the 
pair of parameters that offer the best alignment of the 
object part has to be determined. In order to deal with 
the variances of parameters (a, 6) we let function fs 
act on and on the configuration r = (x, y, a, b). The 
vectors normal to the tangent space of the manifold 
defined by fs{r) = constant are given by ns = 



Thus, the curvature of the manifold tangent space is 
given by the differential 



ns 



dVffs 

IIV.^MI \\S/rfs 



ruldV^fs (14) 



If we project this differential on ns we obtain Ugdns = 
implying that dns lies on the tangent space of 
fs{r) = constant. Moreover, the tangent space of fs{r) 
can be obtained by extending the tangent spaces of 
the projections to (x, y) or (a, b) planes. Namely Is = 

||V^/g|| ( + ^X" ~ ^a^' + ) fs 

where x"" =(0,0, x,?/)^ , x = (x,?/,0,0)^, ^ = 
{y, -X, 0, 0)^, = (0, 0, y, -x)^ = (a, 6, 0, 0)^, a = 
(0,0, a, 6), 13 = (0,0,-6,a)^, (3^ = (-6, a, 0,0). Hence, 
the curvature of the tangent space is expressed via the 

differential vector dns = (^d^,dr^,d^^^ with 

flat representation 

Tr {dns) = (^^(x,a)^^(x,c.)^7 + ^[x,c.]^^[x,«]^^x 

+^l,/3]^^[.,/^]^7f) (15) 

where 7 = ^ + /3, = - c>^^/ 7^ = -^^ + Z^^- 

But the differential vectors dx^ and da^ are related 
with da and dx via their covariance with respect to 
function fsix^a) 



dx^ = — 7: da , da^ = — dx 



da 



dx 



(16) 



where dx-^^ and da^^ are the variances dxr da that are 
normal to the isocontours of fs- 

Covariance of x and a can be obtained by the 
demand that the total differential satisfies dfs = 

d^fsdx + d^fsda = 0, which gives 



da' 



Js 



d.fsdlfs d^fsdlfs ^^^^ 



\\dxfs\ 



dx 



\\dafs\ 



By substituting ( p7| into ([T6| we obtain 



dx" = 



'-^''^^^'-da.da^= 



\\dxfs\ 



Wdafs 



dx (18) 



Equivalently, d^^ and dp^ are related with dp and 
by an analogous concept of covariance between d^ 
and dp with respect to function fs, thus obtaining 



dpfsdjfs 
" ll^xMl' 



dp , dp^ 



ddsdjfs 

" Wdafsf 



di (19) 



Having determined all quantities related with dfis, 
we can return to JTsl) and reformulate it as a dif- 



ferential form using basis df = (drj^dj^dr^^dj^) = 

dr]ns + djli^^^a) + <^^x^[x,a] + ^^^h^^f^]' ^^ich defines 
the volume form dQ = dr] A dj A dr^ A dj^, the line 
form dr = dr]^dj^dr^^dj^ and its Hodge-star dual 
area form dS = ^dT = dj A dr^ A dj^ -dr]A dr^ A dj^ + 
df] Adj A d^^ — df] Adj A dr^ . Next, we define the flat 
curvature of the tangent space as 



where k.^ = J^^^a)^kx,c^)^ = ^^[x,a]' ^ 

Then nsdO. = dfis A ^dV = d{^r]n^df) = 

d{r]djAdr^Adj^), where (i'denotes the exterior deriva- 
tive. Consequently, the integral of the flat curvature 
Ks over a domain by exploitation of the Stokes 
theorem, reads 

/ Ksdn= [ 7?(7,r-,7f)^7A^ir- Ad7f (21) 
Jn Jon 

Namely, the integral of the flat curvature over a 
sub-domain ft in 3?^ is equal to the integral of the 
distances between points of its boundary dft along 
the normals of fs- This property will be used later in 
order to determine the fitting error between a family 
of prototype curves and the image of an object part 
in the curves' family manifold. 

In order to avoid integration of the flat curvature 
over the whole domain we also adopt the inte- 
gral of K'gdVt = dfis A ^dfis = d{^KsngdT), where 
evidently Ks - 



(20) 



Consequently ap- 

plication of the Stokes theorem to the integral of Kg 
over a domain Q gives 

/ K^sdn= [ /4 + + f^\d-f A dr^ A dj^ (22) 
Jn JdQ V 



Consider now the domain ^7 c so as its bound- 
ary is determined by the following three elements: 
1) the interpolated data points' contour rp(^, a, b) 



2) 
3) 



4.1 



of 



(xp(^), yp{^), a, b) as it was obtained in Sect, 
a segment rs{^,a,b) = {xsiO^VsiO^^^^^) 
prototype curve, 

a strip along the curve of parameters 

a{x,y,Sa,P) = (^s{P,x) + 
Sa e [-M«(/3),M«(/3)], p e [po^Po 
where as{P,x) ^ • fsix^c^) = 0}. 
Then, at a point fsiO of the prototype curves sub- 
domain, let the minimal distance of the aligned data 
points sub-domain be (5(^,a,6) = 7^(^,a,6) n5'(^,a,6). 
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The integral of these distances over the prototype 
curves family sub-domain dQs can be obtained using 
pT) and by means of the following formula 



7/(^, a, b)dj A dr^ A dj^ 

2 ' 



(23) 



The problem of finding the optimal fitting of an object 
part data points to a proper parametrized implicit 
prototype curve is now expressed via the demand of 
finding the dfl that minimizes Ei. But, in Appendix 
[a1 it is shown that minimization of Ei demands 
minimization of the integral of equation (|22|, which 
we will call E2{dft). These integrals consider domain 
Q with respect to the manifold defined by the level- 
sets of fsiXi^)- So, under consideration of the basis 
dT together with the relations ( 18 > and {19) the volume 
form of dft reads 



dj A dr^ A dj^ = 



where = \\djsl 




d^ A dSa A dp 



Ad^ Adp 



(24) 



Wdafsl 



doJs 



, d5a 



— n^da and d5^ 



In order to re-evaluate and simplify integration on 
90 using basis (dS^^d^^dSa^dP) we should examine 
how integration on the level-sets of fs{Xi ^) affects in- 
tegration along elements of n^, dS^ and d5a- Namely, 
let equation fs{x^oi) = fc induce two curves, one 
on the X = i^^v) sub-space Xs{S,j(^) and one on the 
a = {a^b) sub-space as{P^x)- Then, in Appendix [b] it 
is shown that if we define the strip along the curve 

[-M,(/3),M,(/3)], 

so as to have constant width, i.e. Ma{(3) = Mq,, V/3, 
alignment errors El{dQ), E2{dft) are simultaneously 
minimized by a pair (^q , /3q ) that satisfies 



\e^{^,p,)\d^dp, = (25) 





T 




ho 


-J- a 





for some e [CS-Co 
(eo,/3o) 



arg mm 

(^o,/3o) 



Lp], and 

(3Q-\-La p^o-\-Lp 



/3o 



^0 



where £v is defined via the formula 



c(C,/?c) = 



Ksd 



a1 



T2 

cx 

X 



o(^-^o,as(^c)) 



rs(C,c.5(/3c)) 



(26) 



(27) 



4.4 Simultaneous determination of the primary 
parameters and tlie segment of the prototype curve 
that optimally fits an object part 

In this section, the results of Sect. |4.3| are exploited 
so as to determine the optimal placement between 









\ Xs(^,0l(N|5)) 


Xs(4,a(P+dp)) \\ 









Fig. 2. Differential domain between an object part 
and a family of prototype curves. The fitting process 
starts by letting an arbitrary prototype curve pass from 
a different point, of the object part each time. In the 
figure the considered object part is the blue line and the 
prototype curve is the red one. Differential movement 
dp along primary parameters curve causes differential 
deformation of the prototype's curvature depicted with 
the differential arc of length d^^ . 



the data points of an object part and the elements 
of a prototype curves family. So, let, once again, the 
considered object part consist of points rp(^i), i = 
1, A/p. In order to evaluate the integrals in equations 
pS] ), ( |26l ), we, first, determine the curve differentials 
dxs and dots along the x and a level sets of fs 
respectively. Namely, dxs = %^<i6 dys = — 
das = -^dp, dbs = ^dp, where, as in Sect. 
4.3 and dp denote the curve length differentials 
along the x and a - level-sets of js respectively. 
Then, starting from the points ao = (<^o, ^o)/ ^§(^0, <^o) 
and for a constant step dp and a given parame- 
ters curve length Lq, we iteratively compute quantity 



4.2 



and by formula 



^xK^n.Pm) using its definition in ( [27| with arguments 

and OLm^x = OLm ^ (^^^'^^^)lrs(^o,a^),a^- Usiugsiug 
Algorithm [l] among all (^o,<^o) we spot the pair 
{xI^olD that satisfies both ^ and (|26|. The fitting 
error between the object part and the prototype curve 
part determined by this procedure is calculated via 
^2(^0^0^^) def ined in S ect. 

ETT{S,T,) = i-^^h{C,.0- 

Among all the examined classes of prototype curves 
and the corresponding stencils obtained by the fitting 
procedure, we decide that each such stencil could had 
generated the object part in hand if the fitting error 
Err{Sjfp) is smaller than a plausible threshold, say 
0.7mm/pixel. In practice, each object part matched 
only one of the potential prototypes with such a small 
error. However, we set an extra criterion for safety, 
namely that if an object part matches two different 
prototypes with such a small error, then we assume 
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that this might have been generated by the stencil 
that offers the smallest matching error. In this way, 
when the above process has been applied to all object 
parts and all potential prototypes, then a one-to-one 
correspondence has been established between any 
object part and a specific potential prototype part. 
However, the exact value of the primary parameters of 
the few stencils that might had been used for drawing 
the wall-paintings, must yet be estimated. 

Algorithm 1 Fitting an object part to a class of proto- 
type curves 

E2 ^ maximal number 



for all ^0 do 
if E2 < then 



c{.in,Pm)\d£,ndp 



end if 
end for 
for all Cc 

Jo 



- E2 

Co 



e [Q,Q + Lj,] do 

m=0 



if = ^0 then 

else if Jo J* < then 

break 
end if 
end for 

return (xS^Q^*) 



4.5 Clustering of the prototype curves that opti- 
mally fit the ensemble of object parts 

Let all extracted object parts with contour point vec- 
tors fp{si), < Si < Lp, p = l,...,Mop as de- 
scribed in Sect. |2.1| and a subset of them Qs attributed 
to the family of prototype curves denoted with S 
via the process described in Sect. 4.4 Moreover, for 
each element rp{si) of Qs^ let the estimated primary 
parameters of the attributed prototype curve be a*. 
We put together all primary parameters estimations 
for all object parts oi Qs forming a set of possible 
parameters vectors for the prototype S over the whole 
wall-painting, SA = [jpeQs 

In order to decide how many stencils correspond 
to the family of prototype curves we classify the 
vectors of SA via an efficient ''k-Means'' method 
[Ts) , thus obtaining Nc clusters of parameters vec- 
tors SAc with minimal Euclidean distances between 
cluster elements. We should also emphasize that each 
cluster SAc of the parameters vectors has a uniquely 
corresponding cluster of the object parts in Q^. 

Now, each cluster SAc corresponds to a poten- 
tial stencil obtained from the S family of prototype 



curves. In order to determine the parameters of these 
stencils and the exact part of them that probably 
guided the drawing of each object part we adopt the 
following process. 

4.6 Exact determination of the primary parame- 
ters of the prototype guides 

In this section we will try to determine the most 
probable pair of primary parameters of a prototype 
that belongs to the class of curves fs{x^y^a^^b^) = 0, 
which is supposed to have been used for drawing 
the object parts that belong to the cluster Q^. In 
order to achieve this, we will use the analysis of Sect. 
4.3 and Appendix [B] for fitting an object part to a 
family of prototype curves, but, this time, with the 
demand that and /3c, which correspond solutions 
of equations (|25| and (|26|, are kept fixed and equal 
to ^p and /3* for each object part p. These values 
correspond to the results of the fitting of each object 
part to the considered class of prototype curves. Thus, 
although ^* varies between different object parts, /3* 
is attributed to zero for obtaining any a* (see Sect. 
|4.4[ ). Thus,we demand that the prototype curve and 
the object part lie on the same primary parameters 
level-set, with initial point a* = 0^5(0, rp(^*)) and 
the prototype curve, Xs{^p,<^p) = rs(^p)/ estimated 
in Sect. HJ 

Having clustered the stencils estimated for each 
object part into groups of prototypes of the same kind 
and of neighboring primary parameters, in this sec- 
tion we look for the prototype curve along which all 
object parts, attributed to stencils of the same cluster 
Qg, can be optimally fitted. Hence, the aforemen- 
tioned demand for fixed primary parameters level-set 
for each pair of a prototype and an object part is now 
violated by asking for a fixed prototype representing 
all object parts of the same cluster. Following the 
analysis of Sect. |4.3| and the definition of domain 
ft between the object part and the prototype, we 
will redefine the strip along the primary parameters 
level-sets a = Q^s'(/3,x)+^a^a foi" given variations 
^ [i^a,Ma] aloug an a - level-set of reference and 
with (3c{Sa) ^ [0, Lq,((5q,)]. Consequently, we redefine 
the domain bounded by the prototype and an ob- 
ject part on the basis of the distances between the 
corresponding curves, S^. Namely, each object part p 
has been already attributed to a stencil part starting 
from point along a prototype curve with primary 
parameters a* and the point ^J(^o) been spotted, 
at which the prototype curve and the object part 
are supposed to intersect, so as the fitting error is 
minimized. But re-positioning the object part along 
a fixed prototype, does not maintain the property of 
^* to be an intersection of the two curves; hence, 
we should re-estimate ^* in [^q^^o + 
correspond to an intersection. Moreover, the starting 
point ^0 of the prototype curve part along which 
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^0 (^c(Q) 




y (pixels) 



index of selected CcC^q) 



Fig. 3. Fitting an object part to a prototype curve. 

Fig.[3a]The optimal fitting error for each tested prototype curve part is depicted with the blue 3D-curve, whose x, 
y coordinates correspond to the starting point of each prototype curve part and its z coordinates correspond to 
the fitting error value. The red planar curve consists of the initial points of each prototype curve part. The point 
that corresponds to the initial point of the prototype curve part that offers the minimal fitting error for the specific 
object part, is depicted with an asterisk. 

Fig. [3b] Fitting error variation Jofe(^^o)) computed via relation (|25| for all relative placements of the considered 
object part and the hyperbola prototypes. Each curve of the figure corresponds to a given initial point along 
the prototype curves, and depicts the evaluation of Jo at each point of the object part. The zero-crossings, 
C(^o)' of these curves correspond to stationary points of the error function, E2. The prototype curve parts (i.e. 
the ^0 selections) that offer evaluations of Jo with a zero-crossing along them correspond to the curves, which 
are colored red. 



^ a 



(29) 
(30) 

In order to evaluate the above minimizer, the points 



arg mm 

^0 



in / / 



a* are ordered along a curve from to aX^, where 
m = argmiiij min^^j \\a* 
Wat 



a* 1 1 , M = arg maxj max^^ ^ 



the alignment error is minimized should, also, be ^c{^x^ that satisfy both to relations below 
re-estimated. Hence, for all possible e [m^,M^], 
where min^ - 

max^{(r5(e,a)-r-;,(0r%(0)}. ^ ^ Ko^^+^p]. 
we define the function ^{S^) which describes each 
possible repositioning of the object part along the pro- 
totype curve's normals so as the domain between the 
prototype curve and the fitted object part is covered 
(see Fig. |4]). 

Using this definition of domain we reformulate 

fin terms of the basis {dS^,d^,dSa,d(3) as in Sect, 
to obtain the fitting error 

E^{p) = / \e^\d^dMS^ 

+ / / / \e''jdpdi,d5'^ (28) 

where functions and are given via 
and Ea respectively, which have been defined 
with equation ( [32| ) of Appendix |B] Boundaries 
determined as domain Q de- 
fines, thus reading = rp{^ — ^0)/ ^ = 
rsi^.aiPc.S^)) and = (x(^c) + ce(/3c, 

In Appendix [C] it is shown that minimization of 
primary parameters estimation error pS) is offered 
by the repositioning ^o(p)/ with intersection point at 



\e^{^,0,5^)\d^d5^ = 



^*||. Then point a* is met at length 5a{p) of 
curve da*{^o{p),Sa) = na{^o{p),Sa)dSa- The proce- 
dure that performs this minimization and returns the 
primary parameters pair {a^,b^) = 0^5(0, r^(^^)) of 
the prototype that optimally fits the object parts of 
is accomplished by Algorithm [2] But, this time, 
in order to obtain an accurate measure for the fitting 
error between each object part and the determined 
stencil part, we calculate the exact values of the point 
by point distances between the corresponding curves 
by means of the lemma introduced in ||4|. In fact, 
for each object part, we properly rotate and translate 
the corresponded stencil part so as their point by 
point distances are minimized, as the aforementioned 
lemma describes. For each such pair of object and 
stencil parts we keep record of the mean and the 
maximal value of their point by point distances. Such 
results for the considered wall-paintings are noted 
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\. \^objectpart 




prototypeV \. 



Fig. 4. Repositioning of an object part along a fixed 
prototype curve. The fitting process runs over all points 
of the object part (in blue) and lets them be the inter- 
section of the object part and the prototype curve part 
(in red). Curve ^{S^) = ^{S^) is eventually selected. 



down in tables Isl [61 



Algorithm 2 Determination of each Object Parts group 
optimal prototype 



maximal number 



for all p e do 

E2 ^ Ea;eQ- Enio \e^{^nAa;)\d^ndS^{p) 
if E2 < e\ then 

end if 
end for 

for all p eQ^ do 

rs{C{p).(^p) '^rs{i{p),a;.) 
end for 

for all ^^GU,,QcKo*(p),a(p) 
for all p e Qs do 



y 



end for 

if p = 1 then 

J*^ Jo 
else if J^r <^ then 

break 
end if 
end for 

return a5(0,rp(^^)) 



n=0 l^x 



Lp] do 

JpdS^ip) 



5 Aspects Concerning the General 
Applicability of the Methodology 

The methodology introduced in this work has been 
motivated by the demands of developing a strict and 
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Fig. 5. [5a) Clustering results for the Hyperbola primary 
parameters determined for the Object Parts of one 
figure of the "Naked Boys" wall painting. The color of 
each group corresponds to the one attributed to its 
optimal guide in Table|4] 

Fig. [5b] Optimal placement of the drawn object parts 
(in blue) along the prototype curves determined as 
optimal representatives of the corresponding object 
parts groups. The object parts of the figure has been 
extracted from the "Lady of Mycenae" wall painting 
and the determined prototype curves correspond to 
the "Hyperbola 2" and "Spiral 2" prototypes as given 
in Table|3] 



quantitative test of the hypothesis for possible geo- 
metric prototypes used for drawing prehistoric wall- 
paintings. However, the approach and the related 
system deal with the problem of determining and fit- 
ting implicit curve models to contour fragments data 
sets in general. Namely, given a set of contour data 
we initially spot the contour segments of the same 
convexity thus forming the set of Object Parts. Then, 
given a set of parametrized implicit curves families 
these Object Parts are attributed to an exact model 
curve and the proper curve segments that optimally 
fit Object Parts are determined. As referred in the 
Introduction, the main advance of the curve pattern 
analysis and fitting methodology developed here is 
that a) optimal relative placement between prototype 
curve and data points and b) the proper form of 
the prototype curve that optimally fits all contour 
segments of its group are simultaneously offered. This 
is achieved by embedding all possible curvature de- 
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formations and all possible, orientation-free, relative 
placements between prototype curve and Object Parts 
into a single error functional. The introduced fitting 
error bares the same minima with the Euclidean one, 
but exploits the functional version of the curvature in 
order to be invariant under rotation and translation. 

One immediate extension of the methodology de- 
veloped in this work, concerns its application to point 
cloud data that do not indicate a contour form; i.e. no 
sense of tangent direction emerges when one reads 
all given data points successively. In the fitting error 
integrals derived by the analysis of Sects. [4Tj |4.1| and 
4.6 the integration is performed along curve length 
- parametrization of the prototype curve. Hence, the 
contour form of Object Parts' data points is induced 
by the fact that for each prototype curve point there 
is one data point met in the direction normal to 
prototype curve's tangent at the specific place. So, 
given a set of point cloud points one could extract 
the contour forms used as Object Parts for the fitting 
process, by treating the point cloud as a compact 
fitting area whose boundary define the curve that 
should be fitted to the prototype curve. Such an 
approach successfully deals with normally distributed 
point cloud around a contour, because the mean curve 
of the boundary coincides with the mean curve of a 
normally distributed set of internal points. This case 
is presented in figure |6| where the contour points of 
an Object Part are distorted independently in x and y 
directions with Gaussian noise, thus forming a point 
cloud. The point cloud envelope is determined by the 
set of points that are assigned as the most distant ones 
from each point of the point cloud in all directions. 
Application of the fitting methodology presented here 
to such point clouds practically offered the original 
prototype curve part. 




Fig. 6. Construction of tine envelope (blue contour) of a 
point cloud (red points) normally distributed around an 
Object Part (black contour in the center). "Mean curve" 
of the envelope (green points) closely approximates 
the original Object Part contour 

A second, not so immediate, extension of the 
methodology concerns application of the fitting error 



formulation and evaluation to problems of estimating 
similarity between shapes. In figure |7| an example 
is presented manifesting application of the fitting 
methodology to estimate optimal matching between 
the shapes of inscribed Greek ''Omega'' - letter sym- 
bols. To such shape matching problems, Euclidean 
distance transform has been used to obtain an implicit 
representation of each shape and the corresponding 
flat curvature values. But for a complete shape match- 
ing application there is a main issue to be resolved. 
Although the introduced fitting methodology offers 
optimal affine deformation of the curvature of a proto- 
type shapes' family so as to match given contour data 
set, the fitting results evaluation relies on the selection 
of the prototype shapes' family. But in the case of 
shape matching analytic form for the prototype shape 
is not available and brute setting of an arbitrary shape 
as prototype does not necessarily offer absolutely 
optimal shape matching results. Hence, following op- 
timal matching between a pair of shapes, one should 
join the matched shapes to a common implicit flat 
representation independently of the selection of which 
shape is considered as prototype. Such a complete 
extension to arbitrary shape similarity problems falls 
outside the goals of the present paper and it could be 
treated in a future work. 




Fig. 7. Evaluation of the fitting error and the resulted 
relative matching for two pairs of inscribed "Omega" 
letters. The figures present the relative matching posi- 
tion between prototype (blue) and data (red) contours 
on the space of the flat curvature evaluations of the 
prototype (the grayscale background). 

Concerning the applicability of the developed 
methodology in cases where hypothesis for existence 
of prototypes in paintings, or in other application's 
contour data, can not be verified and thus it should be 
tested for rejection, we should distinguish 2 cases: a) 
If the implicit forms of the possible curve prototypes 
are given, then in case of rejection each prototype 
optimally fits the corresponding contour data with 
large fitting error. This happens because, even in 
the case of accidental fit between object parts and 
prototypes, there is no common curve that optimally 
fits all object parts attributed to it (see figure [8|. If such 
a curve exists then it is assumed to be a prototype of 
its object parts and the hypothesis is verified, b) If 
the implicit forms of the possible curve prototypes is 



not given, we check if the painting is made by free 
hand as follows: We consider all Object Parts with 
length greater than a proper threshold (depending 
on the size of the painting) and we apply to them 
the methodology pairwise by letting, each time, the 
longer Object Part of the pair to be the prototype 
curve. In this way we obtain a number of classes of 
object parts, where each class consists of the Object 
Parts that fit the assumed prototype with a small error. 
If the number of these classes is large and /or the 
cardinal number of each class is pretty small, then 
we conclude that the considered contours had been 
drawn by free hand. 




Fig. 8. Rejection results for a wall-painting probably 
drawn by free hand. In the left figure the best fit results 
of drawn spirals to the considered model ones are 
shown. Optimal placement of the exponential spiral 
(in red) to the corresponding object parts (in cyan) 
highlights the high error values (3-6 mm/cm). In the 
right figure fitting of the object parts to their common 
prototype, depicts both optimality of the prototype de- 
termination and high error values of the fitting results. 



6 Application of the Methodology 
TO THE Determination of the Method 
USED FOR Drawing Specific Prehistoric 
Wall-Paintings 

6.1 Some historical elements concerning the ex- 
amined wall-paintings 

The so called Mykenaia (Lady of Mycenae), housed in 
the National Archaeological Museum in Athens (inv. 
no 11670,) is the best preserved wall painting on the 
Greek Mainland and of the highest quality. It comes 
from a Late Helladic III (13th century BC) house, near 
the fortification wall, in the area of the Cult Centre at 
Mycenae. 

The majestic face, with the almond shaped eye and 
the elegant profile, has a pensive expression which 
reveals the solemnity of the moment. It is attributed 
to a seated goddess who accepts the offerings of 
the worshippers. In this case it is a necklace made 
of cornelian beads, which she holds tightly in her 
right hand. Alternatively she could be the best pre- 
served of a processional figure bearing offerings to 
a goddess. She wears a short-sleeved bodice over 
a sheer blouse, which delineates her ample bosom. 
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Fig. 9. Photograph of the Lady of Mycenae wall- 
painting housed in the National Archaeological Mu- 
seum of Greece. 

Her intricate hairstyle and rich jewelry (necklaces and 
bracelets) are striking |16|. 

That a device for rendering curves was used for the 
curve of the nose of the Mykenaia was noted early on 

m 

The wall-painting 'Naked Boys' flS^ initially dec- 
orating the internal murals of the lustral basin of 
the edifice 'Xeste 3' excavated at Akrotiri, of the 
Greek island of Thera. The interpretation proposed 
for the scene of the wall-painting is that it depicts an 
initiation rite, during which at least one of the actors 
will achieve manhood. 

.. .. , . . ^ 




Fig. 1 0. Photographs of two Naked Boys excavated at 
Akrotiri, Thera. 



6.2 Choice of possible classes of potential proto- 
type curves based on historical and archaeological 
grounds 

The idea emerged among the authors that certain 
wall-paintings excavated at Mycenae, Thera and Crete 
might have been drawn by usage of prefabricated 
geometric stencils |^|, |Q. A first step to verify this 
conjecture is the determination of a set of geometric 
shapes, whose conception and construction are not a 
priori prohibitive for the era, from an archaeological 
and historical point of view. For example, the linear 
spiral and the hyperbola can be constructed with 
the use of simple tools, even if this might require a 
considerable amount of novelty for the era. 

Thus, extensive archaeological and historical analy- 
sis led to the conclusion that a set of geometric figures 
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TABLE 2 

Descriptive functional forms of the adopted potential 
families of prototype curves 

Model Curve Functional Form 

exponential fEs{x, y\a, b) = -\- y'^ — o? exp (2b6(x, y)) 

spiral 0{x, y) = arctan | + (1 — sgn{x)) ^ 

involute of fiv(x, y\a) = + — (^1 + 0'^{x, y)^ 

a circle Oe + arctan Oe = arctan ^ + (1 — sgn(x)) ^ 

linear spiral fLs(x, y\hi) = x'^ -\- y'^ — hi'^0'^(x, y) 
0(x,y) = arctan ^ + (1 — sgn(x))^ 

ellipse /^^(a;,j;|a,6)=(|)'+(f)'-l 
parabola fp(x,y\a,b) = (ax + b)'^—y 

hyperbola /^(x, 2;|a, 6) = ) ' - (f ) ' - 1 



that could have been conceived and constructed by 
these civilizations, are the following: (a) Exponential 
spiral (b) The spiral generated by unwrapping a 
thread around a peg, usually called the involute of 
a circle (c) The linear or Archimedes' spiral (d) the 
ellipse (e) the parabola (f) the hyperbola. 

It is well known that spiral shapes appear in vari- 
ous prehistoric civilizations even centuries before the 
prehistoric Aegean civilizations. There are infinitely 
many types of spirals. Among them, the involute of a 
circle can be easily generated in everyday life events, 
while the exponential spiral can be found in various 
cockleshells. Thus, not surprisingly rough approxima- 
tions of these 2 types of spirals are encountered quite 
early in various prehistoric civilizations. On the other 
hand, the linear spiral seemingly does not exist in 
nature. In Classical Ages, the conception of the linear 
spiral is so far attributed to Konon in the 3rd century 
B.C.. Next, in ''On Spirals", Archimedes defines linear 
spiral and gives many fundamental properties and 
related theorems. 

6.3 Determination of the stencils and the stencil 
parts most probably used for drawing "Lady of 
Mycenae" and the "Naked Boys" 

In order to test if there are stencils obtained from 
the class of curves defined in Sect. 16.21 we first state 
the equations of these prototypes. Implicit functional 
forms of the adopted prototypes are given in Table |2] 
Analytical expressions for the flat curvature versions 
Ks{x, y, a, b) that correspond to these prototypes fami- 
lies, are obtained straight from these functional forms. 

We have applied the methodology presented in 
Sect. |4] to the prototype equations adopted in Table |2] 
and for the main figures of the wall-paintings "Lady 
of Mycenae" and "Naked Boys". First, photos of the 
considered wall paintings have been taken ("Lady of 
Mycenae") or scanned ("Naked Boys"). These images 
have been processed so as to extract the outlines of 
the painted figures. Next, using the methodology of 
Sect.|2.1[ objects and object parts of these outlines have 



been determined. The object parts with length greater 
than 1cm are employed so as to determine possible 
stencil parts used for creating these object parts. This 
is achieved by means of the procedure described in 
Sect.f 



4.4 



which offers the optimal prototype curve that 
best fits each available object part. After clustering the 
obtained stencil parts into groups of prototype curves 
with neighboring primary parameters values, as Sect. 
4.5 describes, a unique prototype that optimally fits 
all object parts of the same cluster is determined, 
via the methodology described in Sect. |4.6[ We have 
reached the conclusion that both, "Lady of Mycenae" 
and "Naked Boys" wall paintings have been most 
probably drawn by means of linear spirals and hy- 
perbolae stencils /guides. More specifically, it seems 
highly probable that one prototype linear spiral and 
four prototype hyperbolae had been used as guides 
for the drawing of the "Naked Boys" wall painting, 
while two linear spiral and two hyperbola prototypes 
have been used to draw the main parts of "Lady of 
Mycenae" outlines. Equivalently, all object parts of 
the contours of the main figures appearing in these 
wall-paintings impressively match to the prototypes 
specified in Table |3] for the "Lady of Mycenae" wall 
painting and in Table |4] for the "Naked Boys" wall 
paintings. The stencil parts estimated for each object 
part are recorded in Tables |5] and |6] together with 
their lengths and their mean and maximum fitting 
errors. The visual results of fitting the stencil parts 
to paintings' object parts are presented in Figs. 11 
and [l3j We point out, that for space economy reasons, 
the results concerning only two "Naked Boys" main 
figures are presented; very similar results hold for the 
other main figures, too. 

The average matching error of each object part to a 
corresponding part of the prototype stencil, is always 
less than 0.39mm per pixel, while the maximum error 
is always less than 0.86mm. 

7 Conclusion 

In this paper, we have tackled the problem that 
can be described as follows: Suppose that a set of 
outlines of drawn figures is given. Then, look for 
the minimum number of prototype functional forms 
that optimally fit these outlines, with a particularly 
small fitting error. During this fitting process, the 
outlines are partitioned in a number of contiguous 
parts of the maximum possible average length. To 
accomplish the aforementioned, we have shown that, 
in at most two parameters, we call primary, are 
sufficient to describe all curvature deformations of 
a twice continuously differentiable curve. Next, we 
employ the differentiable manifold formed by the two 
primary parameters and the position vectors of the 
considered functional forms. In this 4-manifold we 
perform exhaustive curvature driven optimization, to 
obtain both the optimal primary parameter values 
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TABLE 3 

Stencil Parameters For The "Lady Of Mycenae" 



Type of stencil 


Color 


Primary Parameters (cm) 


hyperbola 1 


pink 


a=11.583 b=12.6025 


hyperbola 2 


brown 


a=6.7500 b=13.7391 


linear spiral 1 


orange 


k=0.0592 


linear spiral 2 


green 


k=0.3175 



Fig. 1 1 . Optimal fitting of the stencil parts of the two 
linear spiral prototypes to the outlines of the "Lady of 
Mycenae" wall-painting. It seems that the two spirals 
were employed for drawing many of the contours of the 
face and hair details of the figure. 



TABLE 4 

Stencil Parameters For The "Naked Boys" 



Type of stencil 


Color in Figure 3 


Primary Parameters (cm) 


hyperbola 3 


magenta 


a=14.24 b=20.12 


hyperbola 4 


light green 


a=4.11 b=6.29 


hyperbola 5 


blue 


a=7M b=17.63 


hyperbola 6 


cyan 


a=2.09 b=2.52 


linear spiral 3 


red 


k=0.169 



Appendix A 

Minimization of the integral of ab- 
solute VALUES VIA THE INTEGRAL OF 
SQUARED VALUES 



and the optimal relative placement of the prototype 
and drawn curves. We would like to point out that 
this fitting optimization behaves particularly well, 
even when the drawings have suffered serious wear. 
Moreover, the introduced methodology is of a more 
general applicability, for example in the case of fitting 
a prototype in noisy point clouds and in optimally 
matching particularly noisy contours. 



Application of this methodology manifests that in 
the fresco technique which was widely spread in 
the Aegean civilizations before 1000 B.C., geomet- 
ric guides were employed corresponding to linear 
(Archimedes) spirals and hyperbolae. It seems that 
this technique had been applied both in Akrotiri, 
Thera in the 17th century B.C. and in Mycenae in 
14th-13th centuries B.C.. However, the parameters of 
the employed stencils differ, as it is demonstrated for 
the first time here. These results indicate perpetuated 
knowledge of constructing these geometric prototypes 
with impressive precision, as the particularly small 
matching error manifests. In addition, it is remarkable 
that these geometric conceptions appear with such an 
accuracy more than 1000 years before their axiomatic 
formulation and treatment by great mathematicians 
such as Archimedes, Euclid, Menaechmos, Appolo- 
nious, etc. One may perhaps speak for an emotional 
sense of geometry in these Late Bronze Age Civiliza- 
tions, which was formally and axiomatically founded 
in the same region in the Classical Ages. 



In this appendix, for the integrals = \c{uj)\duj 



4.3 



we will show 



and /2(^) = Jq c{uj) duj used in Sect, 
that minimization of I2 results minimization of /i. 

In order to determine how minimization of I2 af- 
fects minimization of Ii we exploit the differentials 

dh = \c{n)\dn, dh = c{n)^dn = \c{n)\dii. The 

final relation, together with the one of dli, implies 
that the stationary points of I2 are stationary points 
of /i, even in the case where c{Q) = and that if 
there is a such that WQ, hl^*) < h{^), only 
then Ii{VL*) < h{9). Namely, if there is an Vt' , for 



which h(Q.') < h{Q.*) and if 



then relation dh = 

/ c{Q)dh{Q) < \c 

violates /2(^*) 

does not hold 



sup 

c{n)\dli offers /2(^0 -l2(^*) < 

M\ 



- < 0, which 

< /2(^0- Hence, relation < 



TABLE 5 

Stencil Parts Fitting Results for the "Lady of Mycenae" 



Stencil 
Part 


Type of 
Stencil 


Mean 
Error (cm) 


Max 

Error(cm) 


Max 

Length(cm) 


Hl-11, 










15-19,23 


hyperbola 1 


0.0197 


0.0555 


16.6597 


H12-14, 










20-22,24,25 


hyperbola 2 


0.0169 


0.0449 


9.9723 


Sl-3,5-10, 










18,19,31-34 


spiral 1 


0.0107 


0.0272 


5.2204 


S4,ll-17, 








20-30 


spiral 2 


0.0161 


0.0374 


9.1799 
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Fig. 12. Demonstration of the excellent way, the geometric stencils described in Table |3] fit the corresponding 
object parts for the "Lady of Mycenae" wall-painting. Each color corresponds to a different stencil, as mentioned 
in Table|3] and each determined stencil part is named in correspondence with Table[5| 



TABLE 6 

Stencil Parts Fitting Results for the "Naked Boys" 



Stencil 
Part 



Type of 
Stencil 



Mean Max Max 

Error(cm) Error(cm) Length(c] 



Hb3-5,8,9,11, 

12,14,16,18,23 

Hcl-4,8,9,11,12 

Hbl,2,6, 

Hcl0,13,14 

Hb7,10,13, 

15,20,21 

Hc5-7,15,16 

Hbl7,19,22 

Hcl7,18 

Sbl-8,Scl-20 



hyperbola 3 
hyperbola 5 



0.0209 



0.0168 



hyperbola 6 0.0224 



hyperbola 4 
spiral 3 



0.0193 
0.0295 



0.0559 
0.0477 

0.0470 

0.0555 
0.0497 



14.1541 
6.1036 

7.7096 

7.1760 
15.5649 



X = {x, y) sub-space Xsi.^, and one on the a = (a, h) 
sub-space as{P,x)- Then Xs{i,OLs{P,x)) = Xr and 
€^(3,xs{^,(^)) = (^r implying that ii x ^ Xs then 
dpa = 4^ a = as{Pc) + ^a^a{Pc)r Pc is constant and 
4f-^ e as then d^x = ^ ^ X = Xs{Q + ^x^x(^c), 

is constant. Substituting the volume form |24] of 
in the integral ( [22| ) and separating integrals along xs 
and as, E2 reads 



/3o+I^« M«(/3e) 



Appendix B 

Minimization of E2{dQ) independently of 

PRIMARY PARAMETERS CURVE, as, DEFOR- 
MATIONS 

Let the descriptive equation of S prototype curves 
family fs{Xi^) = fc induce two curves, one on the 



(3o -M^{i3,) ^0 



III K{ic.h.P)\dpd5^di, 



/3o 



(31) 
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Fig. 13. Manifestation of the excellent way, the geometric stencils described in Table [4] approximate the 
corresponding object parts for the "Naked Boys" wall-painting. Each color corresponds to a different stencil, 
as mentioned in Table [4j and each determined stencil part is named in correspondence with Table [e] 



where and £a are defined via the formulas 



Kscf ( 1 + p 




(32) 



and e [^0,^0 + (3^ e [/3o,/3o + Sp^{Q = 

rp{^-^o.a{p,,S^)), r- = (xfe, (5^), -M,(/3))), 
^ = (x(^c,(^x)'^(/^'^a(/^))). If curves rs and 
optimally fit in the sense that their point by point 
distances are minimal, due to continuity of these 
curves there is at least one ^* G [Co, Co + Lp] of zero 
distance between them, i.e. = 0- Then, formula 



pT) implies that there is a point ^^(C), VsiO) of 

zero variance Sa{Pc)- Moreover, consider all possible 
alignments between fs and fp that correspond to 



selecting any point of Cc ^ [Co, Co + Lp] and /3c G 
[/3o,/3o + La], which satisfy C* (3* respectively 
Then, the variance of E2 at each one of these align- 
ments is expressed as SE2 = n . The 

demand of zero variance of E2 at its stationary points 
leads to equations 



^0 l^x 



(C,/3c,o)|dC6^/3c 



(b) 



' rPo- 



\ea{^c.0,/3)\d/3di 



(b)//f:^^"/|:^^n^x(e,/^c,o)idc^/3. 



£«(Cc,0,/3)|dMc 



If we consider the special case where the curve 
Ma{P) that bounds 5^ is a dilated version of as{P^ x), 
then dcx is constant and, because of is equal to zero 
for all p. In this case ea{ic-,P) = and the stationary 
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points (^c,/^) of E2 are given by 



/3o 



\e^{^,Pc,0)\d^dp, = 



(33) 

Hence, expression ( [31) for £^2 is reformulated to 

and hence minimization of E2 is achieved by mini- 
mization of 6E2. Then (Q^Pq) that minimize SE2 read 



iio^Po) = arg min / 

Uo,/3o) J (So 



^0 



(34) 



Appendix C 

Minimization of E2{dn) over primary pa- 
rameters CURVE, as, DEFORMATIONS 

As we mentioned in Sect. |4.6| we look for a 
proper repositioning along the prototype curve's 
normals, where the object part and the prototype 
intersect, thus keeping ^d^^) fixed, or equivalently 
^o(^x) = (^0 + Lp){S^) = ^c{S^). Looking at the fitting 
error ( [28| ), this re-estimation of the intersection points 
induces a selection so as the level Pc of the 
curvature deformation caused by the change of the 
parameters can be kept fixed, say equal to 0. This 
also corresponds to a selection of the proper a* G Qg. 
Among all these pairs of and S^, we seek for 
those that zero the variation of the fitting error 

6E^{p) = . The demand SE^{p) = 

gives rise to the following equations 



(b) /r: ' 



e'^{5^,l3,0)\dm 



'(^,0,5^)|dC< = 



/So 



;(<5^,/3,0)|d/3d<5^ 



where and has been obtained by equation {31) 



with boundaries defined in Sect. 14.61 

But, as we have previously mentioned, the level 
/3 of the prototype curvature deformation along the 
primary parameters level sets is kept fixed and equal 
to and hence integral of vanishes letting the 
E2{p) stationary points equation be 



L 



T 



(35) 



Then expression (28 1 for E2{p) is reformulated to 



5E^{p)dl3, = dp, 



SE2{p) reads 



Q{p) = arg min 



'i^,0,S^)\d^d5^ 



and hence minimization of E2 {p) is achieved by 
minimization of 5E2{p)- Then ^^{p) that minimizes 



^o-\-Lp 



(36) 



^0 
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